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Abstract 

We introduce a model based on a system of coupled nonlinear Schrodinger (NLS) 
equations with opposite signs in front of the kinetic and gradient terms in the two 
equations. It also includes time-dependent nonlinearity coefficients and a parabolic 
expulsive potential. By means of a gauge transformation, we demonstrate that, with 
a special choice of the time dependence of the trap, the system gives rise to persistent 
solitons. Exact single- and two-soliton analytical solutions and their stability are 
corroborated by numerical simulations. In particular, the exact solutions exhibit 
inelastic collisions between solitons. 
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1 Introduction 


The investigation of multicomponent solitons, which arise due to the interplay 
between the second-order dispersion and cubic nonlinearity, has attracted a 
great deal of attention, starting from the classical paper of Manakov [1], and 
further enhanced by works on the copropagation of bimodal waves in nonlinear 
optics 0-0 The dynamics of multicomponent solitons is described by sys¬ 
tems of coupled nonlinear Schrodinger (NLS) equations [see, e.g., recent works 
[8] -[To], and very recent ones mi-lizi dealing with two-component solitons in 
spin-orbit-coupled Bose-Einstein condensates (BECs)]. In particular, the con¬ 
cept of energy sharing in the Manakov model na, or in the modified version 
of this model [H] governed by coupled NLS equations [20] , was a catalyst for 
looking for new integrable models in nonlinear optics [H] , BECs [221I2H] , meta¬ 
materials [SUES], etc. It was found that, in all available integrable systems of 
two coupled NLS-type equations, the ratio between the self-phase-modulation 
(SPM) and cross-phase-modulation (XPM) coefficients, which account for the 
interaction of the components with themselves or with each other, are equal, 
while physically realistic systems depart from this constraint. Therefore, the 
quest for new solvable models involving two coupled NLS equations continues. 
In this context. Park and Shin [26] have developed new forms of integrable 
NLS-type equations going beyond the framework of the conventional Manakov 
model, by adding four-wave mixing (FWM) terms to it. 

In this paper, we investigate a system of coupled NLS equations including 
EWM terms and a time-dependent parabolic potential - generally, with an 
anti-trapping sign, i.e., a potential barrier. An unusual ingredient of the model 
is that signs in front of the kinetic and gradient terms are opposite in the two 
equations. Hence, it does not directly apply to known physical systems. Never¬ 
theless, it is interesting as a “non-standard” nonlinear-wave model. We employ 
the gauge-transformation approach 1271 to construct bright-soliton solutions 
of this system. We conclude that, for a special choice of the trap, bright soli¬ 
tons persist indefinitely long in the system. We verify the analytical results 
by comparing them to the corresponding numerical simulations, and conclude 
that the addition of small perturbations, in the form of sudden variation of 
the trap’s strength, does not destroy the solitons. 


2 The model 


Waves copropagating in optical media interact through the XPM nonlinearity 
[7]. Accordingly, the propagation is governed by the Manakov’s model [T| or 
its generalization 0 - 0 : 
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( 1 ) 


'^Qit + 2^ixx + 2(pii|gi|^ + gi2\Q2\'^)Qi = 0, 

+ 2^2 xx + 2((yf2lklp + fi'22|'?2p)?2 = 0, 

where qj{x,t) [j = 1,2) are envelopes of the field components, and ^22 
account for the strengths of the SPM, while gi 2 and g 2 i represent the XPM. 
It is known that eqs. Q are integrable if either (i) gn = gi 2 = g 2 i = g 22 or 
(ii) 5^11 = 5^21 = —gi 2 = —g 22 - The former choice corresponds to the Manakov 
model proper which has been studied in full detail [TT^dOlldlj . The 

latter choice corresponds to the modified Manakov model [H], in which the 
soliton dynamics has been explored too. 

In addition to the XPM, models of the bimodal light propagation in nonlinear 
birefringent optical fibers include the FWM terms, qlq^ and q^ql, which ac¬ 
count for the coherent nonlinear interaction between two linear polarizations 
of the electromagnetic waves |2l7j . Taking this into regard, we here address a 
novel system of coupled NLS-type equations, including the SPM, XPM, and 
FWM terms with a time-dependent coefficient, and a time-dependent anti¬ 
trapping parabolic potential. The equations are written in the notation corre¬ 
sponding to BEC models based on Gross-Pitaevskii (GP) equations [22123]: 


iQit + ^qixx + l{t){\qi\^ - 2\q2\^)qi --f{t)qlql + ^\'^{t)x^qi^0, 

iq2t + ^q2xx + l{t){2\qi\‘^ - \q2\'^)q2 + l{t)qlq2 + ^A^(t)a:^g2 = 0, (2) 

where 7 (t) is the strength of the FWM terms, and A^(t) is the strength of the 
anti-trapping (expulsive) potential. These potentials occur in various physical 
contexts, such as the interaction of optical and matter-wave solitons with 
barriers |32|-[35|, and splitting of wave packets in interferometers [36] . 

Equations i§ can be derived from the Lagrangian, 


^dqi dql ^ 

^ o ( Oj. Oj. ) o 


dt 


dt 


dqi 


dx 
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1 


+-.mqiV +-.mAq 
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27(^)ki| 1^21 - 


qliql? + qM) 


2/^*\2 


q2 


.dq2 dql 


dt 


q2 


dt 


1 

+ 2 


dq2 


dx 


+ o7(i)k2|^ - -\^{t)x^\q2?, 


with * standing for the complex conjugate. An obvious peculiarity of the 
Lagrangian is that it is sign-indefinite, as the kinetic and gradient terms, 
which contain the t- and ^^-derivatives, respectively, feature opposite signs for 
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components qi and q 2 . For this reason, this system, with “opposite directions” 
of time in the two subsystems, does not describe currently known physical 
settings, although it is somewhat similar to the highly idealized model of 
an optical coupler built of normal and negative-refractive-index cores |37j . 
Nevertheless, the system seems quite interesting in its own right, as a “non¬ 
standard” nonlinear-wave model. 

Another noteworthy consequence of the opposite “time directions” in the two 
subsystems is nonconservation of the usually dehned total norm, 

N=l (\qi{x)\^+ \q 2 {x)\^)dx. (3) 

Indeed, a straightforward corollary of eq. ([^ is the following evolution equation 
for the norm: 

— = 47(t) lm{(ql(x)fql{x)}dx. (4) 

On the other hand, it is easy to check that the system conserves the difference 
between the norms of the two subsystems: 

^/oo “ l^2(a;)|^) dx = 0, (5) 

which is the manifestation of the conservative character of the system with the 
“opposite time directions”. In fact, eq. represents the conservation of energy 
of the dynamical system described by eq.([^. In addition, one can construct 
several conserved quantities as in [38] consolidating the integrability of eq. Q . 
Thus, on the contrary to the “normal” systems, where coherent nonlinear 
coupling leads to exchange of the norm between the subsystems with the 
conservation of the total norm, here the opposite time directions allow the 
coherent coupling to generate or absorb the norm. 


3 The Lax pair 


Equations i§ admits the following Lax-pair representation: 


( 6 ) 
(7) 

where a three-component dost function is $ = ( 01 , 02 , 03 )^, and 


+ = 0 , 

+ = 0 , 
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/ 


u = 


Kit) Pix,t)Kit) Qix,t)Kit)"' 
—Rl{xK)Kit) —Kit) 0 

—R2ix,t)Kit) 0 —Kit) j 


( 8 ) 


/ 


V = 


V21 V22 V23 

X^o \/oo 


(9) 


with 


Vii = —+ iiiit)x(^it) + ^'yit)Ait)Pix, t)Rlix, t)Kit)‘^ 
+ ipit)Ait)Qix, t)Rlix, t)KitY 
i/12 = im^ - Cit))Pix, t)Kit) + "-iPix, t)Kit)), 
i/i 3 = im^ - Cit))Qix, t)Kit) + l^iQix, t)Kit)), 

V 21 = -im^ - t)Kit) + ^(m(x, t)Kit)), 

y22 = Kit)‘^ - i^it)xCit) - i-^it)Ait)Pix,t)Rlix,t)Kit)‘^ 

La 

V<,^ = --QixK)R\ixK)Kit)\ 

= -init)x - Cit))R2ix, t)Kit) + |(i?2(x, t)Kit)), 

Vs 2 = -'-Pix,t)R 2 ix,t)Kit)^ 

V33 = Kit)'^ - icit)xCit) - ^'yit)Ait)Qix,t)R2ix,t)Kit)^, 

with 


P(x,t) = e^^^yP\iix,t) 
Qix,t) = e^yyP">q 2 ix,t) 
RlixR) = 

R2ix,t) = eyy^^P^r 2 ix,t) 
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where 


ri{x,t) 

r2{x,t) 

A{t) 


/ bqi{x,t)q2ix,t)* , dibiq2{x,tfqi{x,t)* 

-aqi{x,t) -;;—;—^-h 


biq2{x,t)qiix,ty 


lit) 


2qi{x,t) 
; K{t) = 


2q2{x,t) 

- cq2ix,t)* + 


2gi(x,t)2 

dbqi(x, t)‘^q 2 {x, ty 


2g2(x,t)2 


1 


Ait) 


; 0(x, t) = fl{t)x^/2 


(10) 


where a, 6, c, d, bi and di are arbitrary constants. One can suitably choose 
these parameters to obtain eq.([^ as the compatibility condition for the Lax 
pair defined by Eqs. (§-@, Ut — Vx + [U,V] = 0, while the spectral parameter 
((t) obeys the following equation: 


with 


Cit) = nit)Cit), 

(11) 

x‘^it) = n^it)-n'it), 

(12) 

^(t) = 

(13) 


It should be mentioned that Riccati equation. (12) has already been employed 
to solve GP-type equations |Hn]-|12]- In fact, the identification of the Riccati- 
type equation (12) gives the first signature of complete integrability of eq. 
([^. Equation (12), which determines the parabolic-potential strength, A^(t), 
demonstrates that it is related to the EWM strength, j{t), through the inte- 
grability condition, which can be derived by substituting eq. (13) in eq. (12): 


- 7"(t)7(t) + 2 iiit))^ - m^\t) = 0. (14) 


Thus, the system of coupled GP equations ([^ (or coupled NLS equations 
with a time dependent harmonic trap) is completely integrable for suitable 
choices of A(t) and 7(t), which are consistent with equation (14). Por constant 
A(t) = Cl, eq. (14) yields 7 (t) = 


It is worthy to mention that the integrable version of eq. (|^ can be trans¬ 
formed, by means of substitution 


Qi,2i^,t) = r— Qi,2iX,T) e:sq){-iflit)x'^/2), 

^/W)m 


(15) 


with X = x/lit), T = Tit), dl/dt = 2Vtl, and dT/dt = 1//^, into a system of 
coupled perturbed NLS equations with constant coefficients. 
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i^ + ^ + {\Qi?-2m}Qi- QlQ\ = ^e(t)Ol, 

+ (2|Q. p - m)Q 2 - QlQl = ie(f)Q 2 . (16a) 

where e{t) = (r2(t) + Thus, if we choose Vt{t) = e(t) = 0 

and the above equation reduces to coherently coupled NLS equation investi¬ 
gated recently in Ref. Q by means of Hirota method. 


4 Persistent solitons and collisional dynamics 

4-1 Analytical results 


To generate bright vector solitons of eq. Q, we now consider the vacuum 
solution (g° = = 0), so that the corresponding eigenvalue problem becomes 




(17) 

(18) 


where 




^-iC{t) 0 0 ^ 

0 iCit) 0 


0 0 iC{t) 


V 


(19) 


/ 


. 


r(o) = 


V 


i({ty + i({t)fl{t)x 0 

0 ~ C,{t)iVL{t)x 

0 0 


0 

0 


\ 


iC,(tY — i(4{t)fl{t)x 
Solving this vacuum linear eigenvalue problem, one gets 


( 20 ) 




$( 0 ) ^ 


dt 


0 


0 


\ 


0 

0 


^iC(t)x+i f Q 

Q gi((t)x+i f ((t)^dt 


( 21 ) 


We now gauge transform the vacuum eigenfunction by a transformation 
function g{x,t) to obtain 
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C/<‘>=9(y<“>g-‘+g.9-‘, 

yd) =9)/(»)9-i+9,9-1, 


( 22 ) 

(23) 


We choose transformation function g{x,t) from the solution of the associated 
Riemann problem such that it is meromorphic in the complex ( plane, as 


'^10 0 ^ 


9ix,t]C) = 




C-Ci 

The inverse of matrix g is given by 


0-10 
0 0-1 




^10 0 
0-10 
^0 0 -1 


\ 

1 

/ 


Cl - cr 

c-cr 


P{x,t) \ , 


(24) 


(25) 


where Ci is an arbitrary complex parameter and P is a 3 x 3 projection matrix 
(p2 _ p-j Pg determined. The fact that and do not develop 
singularities around C = Ci C = Ci imposes the following constraints on 
P: 


P, = {1- P)JU^^\Ci)JP - PJU^^\Ci)J{l - P), (26) 

Pt = {l- P)JV^^\C1)JP - PJV^^\Ci)J{l - P), (27) 


where 


J = 


^10 0 
0-10 
^0 0 -1 


\ 


/ 


(28) 


From the above, it is obvious that one can generate projection matrix P{x, t) 
using a vacuum eigenfunction, $*^°)(x,t) as P = J ■ P ■ J, where 


Trace[M(b] ’ 


(29) 


MP) = $(0)(a),t,c*) .$(0)(a),t,Ci)-C (30) 

In the above equation, is a 3 x 3 arbitrary matrix taking the following 
form 

^*(i)g-2ixi e-2<5i/v^ 0 

0 e~‘^P/\/2 j 

8 


= 


(31) 



such that the determinant of becomes zero under condition 


+ 


-(1)|2 


= 1. Thus, choosing Ci = cti + and using eq. (30), matrix can 


be explicitly written as 


= 


^ e ®^\/2 e e 

e^i/^ 0 


where 


(32) 


0i = 2x(3i{t) ~ ji<^i{'t)(3i{t))dt + 2^1, 

Cl = 2xQ;i(t) - 2 j{ai{tf - /3i{tf)dt - 2xi, 


(33) 

(34) 


with {Q;i(t), /3i{t)} = {ttio, /3io} exp (— / Q{t)dt), while hi and Xi are arbitrary 
parameters. 


Now, substituting eqs. (24) and (25) in eq. (22), we obtain 


[/(') = 


-iCit) [/(o) ^ 

iC(t) 0 
0 iCit)) 


0 ^12^3^ 


2i(Ci - Cl) 


-P 12 0 0 

-P 12 0 0 


(35) 


and similarly for Thus, one can write down the one-soliton solution 


as 


[/(i)^[/(o)_2i(Ci-cr)Pi2, 

y(l)_y(0)_2,(^^_^*)p^3. 


(36) 

(37) 


Thus, the explicit form of one soliton solution can be written as 


= 2Ai£S^^;dosech(di)e*(-«i+^("’*)\ (38) 

= 2A2£^^^;dosech(di)e*(-^i+'^(^'*)), (39) 

where a{t),l5{t) are time-dependent scattering lengths and £ 1^2 are coupling 
parameters, 

with (j){x,t) = Vt{t)x^/2, A\ = A 2 = exp [(1/2) / ^(t)dt], while hi and Xi are 
arbitrary parameters, and £^^^£ 2 ^^ are coupling constants, which are subject 
to constraint |£i^^P + |£ 2 ^^P = 1. 
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Fig. 1. (Color online) Decay of the soliton solntion in the time-independent parabolic 
potential for D(t) = —0.02 (or ^{t) = ea;p[0.02t]), = 0.3, aio = 0.1, /3io = 0.3, 

XI = 0.1, = 0.2. 



Fig. 2. (Color online) Growth of the soliton solntion for f2(t) = 0.02 (or 
7 (t) = ea;p[—0.02t]), = 0.3, aio = 0.5, /3io = 0.5, xi = 0.5, (5i = 0.2. 


Figures and show that, in the case of the time-independent parabolic 
potential, one observes either decay or growth of the bright solitons, for a 
suitable choice of the potential’s strength (r2(t) = const). It should be also 
mentioned that the growth and decay of solitons is a characteristic feature 
of variable-coelRcient NLS-type equations. For example, the density of the 
condensates with exponentially varying scattering length in a parabolic trap 
grows or decays mi with time, depending on the sign of potential, while the 
underlying dynamical system is completely integrable and conservative. 


To stabilize the solitons, we now introduce a time dependence of the parabolic 
potential, selecting Vt{t) as shown in Fig. For this case, the density profile 
of the solution, shown in Fig. indicates that one can sustain the shape of 
the bright soliton. Accordingly, we call solutions of the type shown in Fig. |4] 
“persistent bright solitons”. 
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Fig. 3. (Color online) Evolution of the strength of the parabolic potential, 
A^(t) (which may be both positive and negative), given by Eq. (14), for 
j{t) = exp [(2/3) (1 - 



100 100 


Fig. 4. (Color online) A persistent soliton for the same 7 (t) as in Fig. and 
= 0.3, aio = 0.2, jSw = 0.5, xi = 0.5, (5i = 0.2. 

4-2 Numerical verification 


It is possible to confirm the analytical results by numerical solutions of eq. 
([^, produced by means of the split-step Crank-Nicolson method. In Fig. 
we have plotted the persistent bright solitons derived analytically as per eqs. 
(38) and (39) at t = 10 and 20, and the corresponding numerically generated 
density profiles. Thus, Fig.([^ demonstrates exact matching of the analytical 
solutions to their numerical counterparts. 


Since bright solitons exist due to the special choice of the strength of the 
parabolic potential as a function of time (see eqs. (|T2|)-(0) , we have also 
tested the structural stability of the solitons, by suddenly varying the strength 
of the potential (either increasing or decreasing it by 10%), as shown in Figs. 
^ and (0. From figs. © and 0 , we observe that the addition of a small 
perturbation does not impact the stability of persistent solitons. 
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X X 

Fig. 5. (Color online) Comparison of analytical and numerical solutions for solitons. 
(a) The upper left panel: the qi component ai t — 10; the upper right panel: q 2 at 
t = 10. (b) The bottom left panel: at t = 20; the bottom right panel: q 2 aXt — 20. 
Parameters are the same as in Fig. @- 

5 Collisional dynamics of bright vector solitons 


The gauge-transformation approach can be easily extended to generate mul- 
tisoliton solutions [27] . In particular,, the two-soliton solution q^l for the two 
modes can be expressed as 


gf) - 2IAqB, gf = 2 IA 2 /B, (40) 

where 


Ai — M121M222 (C2 — Cl) (Cl ~ Cl) (C2 — C2) + M122M221 (C2 — Cl) (C2 — Cl) (C2 — C2) 
+ MniMi22 (C2 - Cl*) (C2 - Cr) (C2 - C2) + M.uMui (Ci - Ci*) (C2 - Ci) (C2 - CD , 
A2 = M112M2U (C2 “ Cl) (Cl “ Cl) (C2 “ CD + -^111-^212 (C2 “ Cl) (C2 “ Cl) (C2 “ CD 
+ M212M221 (C2 - CD (Cl - CD (C2 - C2) + M211M222 (Cl - CD (C2 - Cl) (C2 - CD, 
B = {M122M211 -\- M121M212) (Cl “ CD (C2 — CD + (-^112-^221 + M111M222) (C2 ~ CD 
(C2* - Cl) + (MniMn2 + M221M222) (C2 - Ci) (C2* “ CD , 
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Fig. 6. (Color online) Density profiles of the fields produced by suddenly increasing 
the potential’s strength for, = exp[(20/3)(l — at (a) t — 10, (b) t = 20. 

Mio, = 


Miij = e ^^V2; ivii2j = a ■■'fci ; iwisj 

M2ij = M22j = e®Vy2; M2Sj = 0; 

62^^ ] M^2j = 0 ; Mssj = 


-2 5 


dj = 2xPj{t) “ 4y {aj{t)l3j(t))dt + 25j, (41) 

Ci = ‘2xaj{t) - 2 j- ^j{tf)dt - 2xi, (42) 

and j = 1, 2 

In Fig. one can observe inelastic collision of persistent solitons. The colli- 
sional dynamics predicted by the analytical solution (the top panel in Fig. 


13 
















100 



X 



-1-1 —I-1 "I I 

- 3 - 2-10 1 2 3 


X 

Fig. 7. (Color online) The density profiles of the fields obtained while suddenly 
decreasing the potential’s strength, for 7 (f) = ea:p[(l/15)(l — at (a)f = 10, 

(b) t = 20 . 

and its numerical counterpart (the bottom panel in Fig. are identical. 


6 Conclusion 


The aim of this work is to investigate the dynamics of solitons in the inte- 
grable system of coupled NLS equations with “opposite directions” of time 
in the two subsystems. The system includes the time-dependent nonlinearity 
coefficient, which must be specifically related to the coefficient in front of the 
parabolic-potential terms, to secure the integrability. By means of the gauge 
transformations, we have demonstrated that a special choice of the time de¬ 
pendence of the trap may effectively stabilize bright solitons. We have also 
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Fig. 8. (Color online) Inelastic collision of solitons for the choice 
7(t) = exp [(2/3) (l-e-O-3*)]., aio = 0.1, ^20 = 0.25, Ao = 0.3, = 0.2, 

= 0.1, ^2 = 0.2, XI = 0.3, X 2 = 0.4, = 0.85i, = O.Ssuch that 

+ |e 2 ^^P = l,(j = 1,2). The top and bottom panels show the analytical 
solution and its numerical counterpart. 

observed inelastic collision of persistent solitons for the same choice of trap 
frequency which is subsequently confirmed by numerical simulations. 
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